Dilute Bose gas: short-range particle correlations and ultraviolet divergence 
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The modified Bogoliubov model where the primordial interaction is replaced by the t matrix is 
reinvestigated. It is shown to provide a negative value of the kinetic energy for a strongly inter- 
acting dilute Bose gas, contrary to the original Bogoliubov model. To clear up the origin of this 
failure, the correct values of the kinetic and interaction energies of a dilute Bose gas are calcu- 
lated. It is demonstrated that both the problem of the negative kinetic energy and the ultraviolet 
divergence, dating back to the well-known paper of Lee, Yang and Huang, is connected with an 
inadequate picture of the short-range boson correlations. These correlations are reconsidered within 
the thermodynamically consistent model proposed earlier by the present authors. Found results are 
in absolute agreement with the data of the Monte-Carlo calculations for the hard-sphere Bose gas. 
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I. INTRODUCTION 

The observations of the Bose-Einstein condensation in 
the magnetically trapped alkali atoms Q significantly re- 
newed interest in the theory of the Bose-Einstein con- 
densation (see, e.g., Ref. ||) and also motivated its ex- 
tensive re-examinations. The standard method, which 
is widely used to study the Bose gas of neutral atoms, 
operates with replacing a pairwise interatomic potential 
by the effective, or "dressed", one (this is why below 
the name "effective-interaction" is used for the approach 
whose various representations are listed in Ref. ||). Such 
a replacement allows one to overcome calculational ob- 
stacles coming from the fact that realistic potentials are, 
as a rule, strongly singular pL For a dilute Bose gas 
this method, in its simplest form, is reduced to the re- 
placement of the real potential by the zero-momentum 
t matrix obtained from the ordinary two-body problem. 
In so doing, a new difficulty arises. Namely, the so-called 
ultraviolet divergence appears in the perturbation series. 
This feature is usually considered as a nonfundamcntal 
one, and textbooks do not draw much attention to it. 
However, some troubles compel us to re-investigate this 
problem in more detail. 

It is well-known that the first microscopical treatment 
of the Bose-Einstein condensation in an interacting Bose 
gas has been realized by Bogoliubov in his classical pa- 
per of 1947 H and is concerned with the weak-coupling 
regime. In the Bogoliubov model, single-particle excita- 
tions coincide with the collective ones, the latter can be 
obtained within the dielectric formalism in the random 
phase approximation, which is reduced to summation of 
infinite number of the bubble diagrams j(],|7j . In order to 
generalize the Bogoliubov model to the case of a strongly 
singular potential one should take account of the multi- 
ple two-particle scattering, and, thus, summarize infinite 
number of the ladder [t matrix] diagrams. The result 
of this summation is considered to be expressed in the 
replacement of the Fourier transform of the pairwise po- 
tential $(fc) by the t matrix t — 4irh 2 a/m, a being the 



scattering length and m denoting the boson mass. How- 
ever, it is the replacement that leads to the ultraviolet 
divergence. This can easily be shown with the help of 
the expansion of the energy per particle for a weakly in- 
teracting Bose gas in powers of the gas parameter na^ 
(n = N/V is the boson density): 
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Here clq and a\ are the leading and next-to-leading terms 
in the Born series for the scattering length a given by 



a = ao + a% + a-i + ■ 
m 
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with Tk — h 2 k 2 /(2m). For more detail see the article of 
Brueckner and Sawada in Ref. || and paper ||. As it 
is seen from Eq. (||), the replacement <!>(fc) — > Airh 2 a/m 
yields ao — > a and a\ — > oo. This problem is consid- 
ered to be solved since the PhD thesis of Nozieres (1957) 
where it has been found that the divergence is an artifact 
of combining the bubble and t matrix diagrams and 
appears due to double account of the term 2irh 2 nai/m (a 
already contains ail). We emphasize that there is no dou- 
ble account in the original Bogoliubov approach (ao does 
not include ax). Thus, replacing original interaction <E>(fc) 
by the t matrix, one should be careful and ignore the di- 
vergent term proportional to J d 3 k/k 2 . This allows for 
obtaining the classical result of Lee and Yang [jj) given 

by 



2irh 2 na 
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and first found in the framework of the binary collision 
expansion method. Note that in the pseudopotential for- 
mulation escaping the divergence corresponds to using 
47rfi 2 a d _.. . . , rTr A?:h 2 a a 

V eff = — — Q^ rS ^ mstead 01 Veff = —^8{r) 0. 
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The present article shows that the trouble that mani- 
fests itself in the form of the ultraviolet divergence, dis- 
cussed above, can not be cured by removing the diver- 
gent term. Indeed, in the next section the kinetic and 
interaction energies of a strong-coupling Bose gas are in- 
vestigated. It is shown that the modified Bogoliubov 
model, where the original pairwise potential is replaced 
by the t matrix, yields incorrect values of these impor- 
tant quantities. In particular, the kinetic energy turned 
out to be negative. Note that the subject of this section 
is of interest not only in the particular context of the 
effective-interaction scheme. To the best of our knowl- 
edge, the interaction and kinetic energies of a strongly 
interacting Bose gas have never been investigated in de- 
tail. Moreover, one can find absolutely different points 
of view concerning these quantities in literature. For in- 
stance, there is opinion that the mean energy of a Bose 
gas, taken in the leading order in the gas parameter na 3 , 
is equal to the interaction one |10|. According to another 
point of view |ll[] , all the mean energy of a Bose gas of 
hard spheres is kinetic in the same order. Now there 
is essential need of clarifying this ambiguous situation 
because it directly influences interpretation of the exper- 
imental data (see the first paper in Ref. [^0|). Further, in 
Sec. [II the pair distribution function g(r) is calculated 
within the effective-interaction approach. This quantity 
turned out to be unphysically negative at small boson 
separations. Remind that g(\ri — Y2\)/V is the density 
of the conditional probability of finding one particle at 



the point ri while another is at the point T2- In Sec. IV 



a correct expression for g(r) is derived with the help of 
the Hellmann-Feynman theorem and a variational the- 
orem for the scattering amplitude. To reveal a reason 
of the failure of the effective-interaction scheme, we con- 
sider a representation for g(r) in terms of the in-medium 
pair wave functions. On the basis of this formula one 
can conclude that the Bogoliubov model modified by the 
replacement $(fc) — * Anh 2 a/m does not fit the strong- 
coupling regime due to weak-coupling features that sur- 
vive in the effective-interaction approach. A possible way 
of the strong-coupling generalization of the Bogoliubov 
model, which is free from the troubles of the effective- 
interaction scheme, is considered in Sees. M, (v| and VII . 
It reproduces the result of Lee and Yang (m and, in addi- 
tion, yields correct results for the kinetic and interaction 
energies and pair distribution function. Note that we deal 
with a systematic calculations in the first two orders of 
the expansion in na 3 . Calculations of the next terms im- 
ply additional and rather extended investigations beyond 
the scope of the present paper. By the moment, the au- 
thors are able to contemplate some important points of 
studying the strong-coupling regime only in terms of the 
in-medium pair wave functions. How to proceed with 
this interesting problem in the framework of the more 
familiar Green's function method is an open question. 



II. THE INTERACTION AND KINETIC 
ENERGIES 

Let us start from the original Bogoliubov model and 
find the expansions of the kinetic and interaction energies 
in powers of the weak-coupling gas parameter nap 5 . The 
most simple way of doing this is based on the well-known 
Hellmann-Feynman theorem given by 



SE = (SH). 



(5) 



Here (• • •) = (0| ■ • • |0) stands for the ground-state aver- 
age, SE and SH are infinitesimal changes of the ground- 
state energy E = (H) and the Hamiltonian 



H 



E 



2m 



(0) 



where 7 is an auxiliary parameter, the coupling constant, 
which should be put equal to unity in final formulae. 
According to this theorem for the kinetic and interaction 
energies per particle we have 



d 3 r 7$(r)g(r) = 7 
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where e = E/N = Ekin + tint is the total energy per par- 
ticle, nk denotes the occupation number. From Eqs. (|l]), 
©> (0) an d © on e can derive for the original Bogoli- 
ubov (B) model the following equations: 
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We remark that a\ < and, therefore, Ekin in Eq. ( |10| ) 
is positive. One can see that Skin -C £mt (|ai| *C ao!) 
and, thus, the main part of the mean energy comes from 
the boson-boson interaction in the weak-coupling reg ime. 
It should be stressed that the formulae d) and (|TJ) can 
be obtained directly from the Bogoliubov expressions for 
Uk and g(r), since the original Bogoliubov model is fully 
consistent. Now what we have for the modified Bogoli- 
ubov model with the original interaction $(fc) replaced 
by t — Airfi^a/m? Substituting a for ao and removing 
the terms depending on 01, one can readily obtain for 
the modified Bogoliubov approach (mB) the following 
relations: 
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It is known that a should be positive, otherwise the Bose 
gas would be unstable. In this case s\™^ in Eq. ([t2Hs 



negative. Thus, the result given by Eqs. ( |11| ) and 
can not be correct. 

What are the true values of Ekin and Si n t in the strong- 
coupling regime? This can again be clarified with the 
help of the Hellmann-Feynman theorem (|5|). However, 
first of all one should prove a useful variational theorem: 



<5t/ (0) (0) 



d 3 r 
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+V {0) {r)6(<P{r))^°\r)~\, (13) 
where we define the scattering amplitude by the relation 



[/ (0) (fc)= / d 3 r^ (0) (r)$(r)exp(ikr), 



(14) 



and ip^ (r) = 1 + (r) is the solution of the two-body 
Schrodinger equation in the centre-of-mass system 



VV 0) (r) + $(rW 0) W = 

TO 

with the asymptotic behaviour at r — ► oo 
(p(°\r) ->l-a/r. 



(15) 



(16) 



From Eqs. (|T|) and (@ it follows that U^(0) = 
Airh 2 a/m. Relation (pL3|) can be proved representing 
Eq. (Q) in the form 



U (0 \0) 



d 3 r 



-|V^ (0) (r)| 2 + [^ (0) (r)] 2 $(r)l, (17) 



using integration by parts, and taking into account the 
Schrodinger equation (|l^) and the boundary condition 
(p"6|). Further, varying Eq. @, we get Eq. @ with 
Eqs. ( |l5|) and (|l6| ) taken into account. 
From the theorem ( |l3] ) it follows that 
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da 



7— = to— = a 
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d 3 r [<^(r)] 2 $(r), 



where we introduce by definition the quantity 

i2 



b =^J d 3 r|W (0) (r)| 



(18) 
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As it is seen, Eqs. (|J 
lowing expressions: 

_ 2-!Th 2 (a~b)n 
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@ and (|T|) provide the fol- 
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Equation (|T^) implies that & is a positive quantity that 
can be considered as a new characteristic length. We 



stress that b is expressed in terms of a and its derivatives 
[see Eq. (|l8|)] and determined by a shape of the inter- 
action potential. For example, when <i>(r) is the hard- 
sphere potential 



$(r) 



+oo, r < a 
0, r > a, 



(22) 



we have b = a, which follows from the obvious relations 
da/dj — da/ dm = 0. While for the weak-coupling po- 
tential Q b ~ — ai, a ~ ao, and, hence, b ^ a. The 
latter implies that ei n t ~ £ and Ekin ~ in the weak- 
coupling regime. This allows for concluding that Ei n t 
and Skin depends on a particular shape of the interac- 
tion potential even in the leading order of the expansion 



in na 3 . Equations ( |20|) and ( [2l|) te stify that in the case 
of the hard-sphere interaction (|22j) all the mean energy 
is kinetic, which agrees with the expectation of Lieb and 
Yngvason |ll|] and definitely excludes the variant adopted 
in Ref . jl0| . One can see that this result is rather general: 
for the hard-sphere potential (2i) the interaction energy 



is equal to zero for any density. Indeed, $(r) given by 
( p2|) can be thought of as the limiting case of the repulsive 
potential 



$(r) = 



Vq, r < a 
0, r > a. 



It is clear that saturation takes place when Vq S> e: fur- 
ther increase of the parameter Vq does not change the en- 
ergy per particle e. Hence, according to Eq. (Q), £ int = 
because ds/d^f — at 7 = 1 in the limit Vq — > +00. 
Thus, the incorrect results of the effective- interaction 
scheme given by Eqs. ( pT| ) and ([l2]) are, say, of the weak- 
coupling character as they imply the relations £i n t ~ £ 
and Ekin ~ in the leading order. 



III. PAIR DISTRIBUTION IN THE 
EFFECTIVE-INTERACTION APPROACH 

It turned out that the trouble concerning the inter- 
action and kinetic energy is accompanied within the 
effective-interaction approach by the problem related to 
the pair spatial correlations. Let us show this with the 
help of the Beliaev's paper and article of Hugenholtz and 
Pines. In the latter the structural factor S{k) has been 
calculated. We can employ it to find the pair distribution 
function g{r) with the well-known relation: 



9(r) 



1 



1 



d J fc 

(27T) 5 



[S(k) - l]exp(tkr). 



(23) 



However, Hugenholtz and Pines ||| restricted themselves 
to the approximation valid at k — > 0, which makes it im- 
possible to integrate their result to get g(r). So, their 
calculations are repeated below but not involving the ar- 
gument of k — > 0. They investigated S(k) with the help 
of the Green function 



3 



F k (t - t') = -^(0|T{p k (i)p_ k (t')}|0), (24) 

where = a^a q +k is the Fourier transform of the 
density operator. According to the paper of Hugenholtz 
and Pines ||, in the lowest order in the condensate de- 
pletion (n — no)/n (when one can take n = no, no is the 
density of the condensate) the quantity F k (t — t') can be 
written as 

F k (t-t') 

= -ijV<0|T{oL k (t>+^k(t)}{4(* , >^-k(* , )}|0), (25) 
so that its Fourier transform is given by 
F(k,w) = 

N [2T k + E+ + E n -2fi- 2E 02 ] 



(E+ - E n )] - [T fc -/i+| (S+ +S n )] + S 



02 



-*<J 

(26) 



where Ej^Sr^ and E02 are the effective potentials in- 
troduced by Beliaev; fj, stands for the chemical potential. 
Beliaev has found that in the lowest order in the con- 
densate depletion the effective potentials and fi obey the 
following relations: 

H = n/(0, 0), E 02 = n/(k, 0), E± = 2n/.(k/2, k/2). 

(27) 

Here the expression /(k, p) is the "non-diagonal" scat- 
tering amplitude: 

/(k,p) = y"d 3 r^°)(r)$(r)ex P Hkr), (28) 

where ipp (r) obeys the Schrodinger equation ( |l5| ) with 

the right-hand side equal to (h 2 p 2 /m)ipp\r) Be- 
sides, / s (k,k') = Jf(k,k') + /(-k,k')]/2. Substituting 
Eq. ( p7j ) in Eq. (^6|), we arrive at 

2N [T k +n(2f s (k/2, k/2) - /(k, 0) - /(0, 0))] 



with 



ft 2 (w — + id) (u> + LUk — id) 



fkv k = ^(T fc +2n/ fl (k/2,k/2)-n/(0,0)) -n 2 / 2 (k,0). 

Now, to derive the structural factor 5(fc), one should 
utilize the definition ( p4| ) at t — t': 

+00 



After transparent integration one can find 

_ ^ + "(2/,(k/2,k/2) - /(k,0) - /(0,0)) 



With the help of the limiting relation at k — > 

47T?7 2 <7 

/.(k/2, k/2) ~ /(k,0) ~ /(0,0) = 



equation (29) is reduced to the familiar expression 

S (k) ~ i*. ~ Tk = fc2 

hw k y/j® + 2nT fc /(0, 0) V/c 4 + 167ma/c 2 

(30) 

As it has been mentioned before, Eq. ( p0| ) can not be 
used in Eq. (|||) because the integral diverges in this sit- 
uation at large momenta. On the contrary, the exact 
(in the lowest order in the depletion) variant ( p9f ) can be 
employed without problems concerning integration. In- 
serting Eq. (p9|) in Eq. (23), at n — > we get 



d 3 fc /(k,0) 
(2tt)3 T k 



exp(ikr). 



(31) 



Equations ( |l5| ) and ( p8| ) at p = allow for representing 
Eq. (|3l| ) in the limit n — > as 



lim c/(r) = l + 2^°\r). 

n—*0 



(32) 



Now we can easily be convinced that the result for g(r) 
derived within the effective-interaction scheme is inade- 
quate in the strong-coupling regime [0]. Indeed, in this 



(0) 



(r = 0) = and, by definition, f/>(°)(r 



regime tp 

0) = —1. This, taken together with Eq. (|32]), leads to 
g(r = 0) — ► — 1 when n — ► 0. Thus, the pair distribution 



function inferred from Eq. (29) is negative at small bo- 
son separations, which implies unphysical picture of the 
short-range boson correlations. 

In addition, one can see from Eq. ( p2| ) that the 
effective-interaction approach is not self-consistent. In- 
deed, in the case of a strongly singular potential we have 
g(r = 0) — ► — 1 at n — ► 0. Then, from the first equality 

in Eq. (fj]) and Eq. ( fj2| ) one can find to be infi- 

nite (one more divergence!). However, Eq. ( pd| ) gives the 

finite value for . Thus, two different ways of calcu- 

lating the interaction energy yield two different results. 



IV. PAIR WAVE FUNCTIONS 

Now the question arises: what behaviour is correct for 
g(r) in a dilute strongly interacting Bose gas? This can 
be clarified with the help of the theorems (|J) and (|l^) , 
starting from the correct expansion (fl|). For a homoge- 
neous system the pair distribution function can be rep- 
resented as 



7(r) 




4 



Here S(r) stands for the Dirac's delta-function. 
Hellman-Feynman theorem (^) yields 

2 Se 
9{r) = -- 



The 



n 5$(r) ' 

where we put (N— l)/V ~ n in the thermodynamic limit. 
Substituting the expression (Q) and using the theorem 
131), we arrive at 



9 (r) = [^°\r)f 



(33) 



This equation is discussed in detail in Sec. [v]. In the 
limit n — ► from Eq. ( |33"1) we obtain the expression 
known since the Bogoliubov's article [see the conclud- 
ing part of his classical paper where a possible way 
of going beyond the Born approximation is discussed] : 



limo(r) 

n— ► [) 



/'(r)] =(l + ^(r))' 



(34) 



Here <p(°\r) obeys the Schrodinger equation (|l5| ) with 
the boundary condition dig ). To compare the results 
given by Eqs. ( |32|) and (|34|), let us consider the simplest 
case of the hard-sphere pairwise interaction defined by 
Eq. (§2|). In this case we get 



tp<®(r) = 



0, r <a 

1 — a/r, r > a. 



(35) 



The data found from Eqs. and @ for the hard- 
sphere bosons are presented in Fig. [I]. As it is seen, at 
small particle separations r < 3a the difference between 
the curves (1) and (2) is essential, while at large r these 
curves are close to one another. We emphasize that the 
incorrect behaviour at the short distances is related to 
the modified Bogoliubov's model [with the replacement 
<£>(&;) — > 2irH 2 a/m] but not to the original one developed 
for the case of a weakly interacting Bose gas. 

The failure of the effective-interaction approach can be 
understood with the help of the interesting relation con- 
necting the pair distribution function with the in-medium 
pair wave functions and following from the Bogoliubov 
principle of the correlation weakening |l3|-|l5|]. We re- 
mind that the pair wave functions, by definition, are the 
eigenfunctions of the reduced density matrix of the sec- 
ond order. For the Bose gas, a system with a small con- 
densate depletion [n — no)/n <C 1, the pair distribution 
function can be written as 

g r = f — ) ip r +2— / jt-t q /2 (r , 36 

which is accurate to the next-to-leading order in (n — 
rio)/n [T^ ]. Another restriction for this representation is 
an assumption that there are no bound pair states (see 
details in Ref. 13|). In Eq. (Bq) <p(r) is the in-medium 
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FIG. 1. The pair distribution function for the hard-sphere 
bosons in the limit n — > 0: 1 - the data of Eq. (jj^l); 2 - the 
result of the effective-interaction approach Eq. (^2|). At 
r < 2a Eq. (Jf^) yields unphysically negative values for g(r). 

the quantity y q /2(r) denotes the in- medium wave func- 
tion of the relative motion of the pair of bosons with the 
total momentum S,q. This pair includes one condensed 
and one uncondensed particle. The functions <p(r) and 
<^ q / 2 (r) are chosen as real quantities. It is also conve- 
nient to introduce the in-medium scattering waves ip{r) 
and fp p (r) (p ^ 0) given by 

(p(r) = 1 + ip(r), <p p (r) = V2cos(pr) + ^ p (r) (37) 
with the boundary conditions at r — ► oo 



ip(r), V P (r) -> 0. 



(38) 



The Fourier transforms of the scattering waves can be 
expressed in terms of the Bose operators aj, and a p [h_3| : 



n 



V 
2n 



\ a 2p a P+k a p-k; 

n 2p 



(39) 



At n — > the in-medium pair wave functions tend to 
solutions of the ordinary two-body problem. For exam- 
ple, <p(r) — > </? < -°- ) (r). Hence, as the condensate depletion 



pair wave function of two condensed bosons; n q — (a q a q ) 



(n — n a )/n approaches zero when n — ► 0, Eq. ( |36[ ) is 
reduced to the Bogoliubov's equation (|34|). A weakly in- 
teracting Bose gas is characterized by a minor role of the 
particle scattering, so that \ip(r)\ <C 1 an d [V'pWI 1- 
In particular, the Bogoliubov model of a weak-coupling 
Bose gas implies ip p (r) = [jl]|[l4|, while tp(r) ^ 0. Thus, 
within the Bogoliubov model, Eq. ( |36| ) can be approxi- 
mated as 



5 



g{r) = 1 + 2ip(r) + 



d J fc 

(27T) 5 



n fc exp(ikr), 



(40) 



where the terms of the order of (n—no) 2 /n 2 and ip(r)(n- 
no)/n have been ignored. At n —> one can find P,H 



. 9 (r)^l + 2v4 0) (r), 
where ipg = ip^\r) — 1 obeys the equation 

- — VV ( B 0) (r) + $(r) =0, 
m 



(41) 



(42) 



which is the two-body Schrodinger equation (|Tq) taken 
in the Born approximation. We stress that Eqs. ([40|)-(|42]) 
are related to the original Bogoliubov model. As to the 
modified Bogoliubov approach with the t matrix, it pro- 
duces Eq. ([32]) . Now we are able to explain the failure of 
the effective-interaction scheme. The replacement of the 
original interaction $(fc) by the t matrix leads to passage 

from tph (r) to ^°\r). However, it does not influence 
the Bogoliubov ansatz for the relation between the 
pair distribution function and the pair wave functions. 
Substituting ip(k) — (aka_k) /no and n k = (a^ak) in 
Eq. ( p0| ) and using Eq. (|2^) , one can see that the approxi- 
mation ( prCf ) is consistent with the approximation ( g5| ) for 
the Green function (p4|). Thus, the effective-interaction 
scheme combines the features of both the strong- and 
weak-coupling regimes, and this is the actual reason for 
its problems discussed in the previous sections. In partic- 
ular, the ultraviolet divergence does come from the dou- 
ble account mentioned in Sec. |, but our analysis reveals 
that the double account, in turn, is a consequence of the 
combination just pointed out. Note that the thermody- 
namic inconsistency of the modified Bogoliubov model, 
found in our previous publications P,p^4|, results from 
this combination, too. 



ordinary two-body problem at n — > 0. The model consid- 
ered below is similar to the Brueckner theory (see Sees. 
36 and 41 in Ref. |l7j) but with one advantageous excep- 
tion. The point is that the Brueckner theory implies the 
momentum distribution of interacting particles to be the 
same as in an ideal gas. As to our model, it results in a 
system of equations connecting the in-medium pair wave 
functions with the momentum distribution. 

In addition to Eq. fl36|), the proposed strong-coupling 
generalization of the Bogoliubov model involves two 
other relations. The first is the familiar expression for 
the mean energy per particle 



d J fc 



Ilk: 



T k — + - d J r ff (r)$(r), 



(43) 



which can be found in any textbook on statistical me- 
chanics. The second is, say, the semi-phenomenological 
relation that connects the scattering waves ([37]) with the 
momentum distribution n k : 



nk(nk + 1) = nlijj 2 (k) + 2n 



(2tt) 3 



n 9 < /2 (k). (44) 



It is worth saying several words about Eq. ([44]). This 
equation represents, in particular, the well-known fact 
that when there is no scattering [interaction] in the sys- 
tem, there are no uncondensed bosons n k = Q. The 
larger interaction, the larger depletion of the Bose con- 
densate. Equation ( [4~4| ) generalizes the similar relation of 
the Bogoliubov model given by (see Refs. Iqjl4||) 



nkijik + 1) = n 2 ip 2 (k). 



(45) 



The generalization (Q) has been chosen in pjli}] for 
the following reasons. First, it provides the relation 
^/2tp{r) — linip^o ?A p (r) as well as 



V2ip(r) = lim y p (r), 



(46) 



V. STRONG-COUPLING GENERALIZATION OF 
THE BOGOLIUBOV MODEL 

It follows from the consideration given above that in 
the strong-coupling regime one should leave the effective- 
interaction approach and develop a new one. What way 
should one prefer? At present the authors have no final 
solution providing the solid theoretical scheme and mak- 
ing it possible to realize systematic calculations like in 
the weak-coupling case. However, we can propose a rea- 
sonable strong-coupling generalization of the Bogoliubov 
model based on Eq. ( p6[ ) and semi-phenomenological re- 
lation (Q) discussed below. This generalization is jus- 
tified as it reproduces the result of Lee and Yang (||), 
gives correct picture of the short-range particle correla- 
tions, and yields Eqs. (20) and (pl[). Besides, the way 
proposed provides self-consistent calculations of the in- 
medium pair wave functions approaching solutions of the 



which can be inferred from Eqs. ( p7| ) and (|38|). Second, 
it leads to the correct thermodynamics and behaviour 
of g(r) at short distances for a dilute Bose gas in the 
leading and next-to-leading orders in no 3 . Third, the 
relation ( ^i| ) results in Eq. ( p4| ) that gives not only short- 
range but also correct long-range behaviour of f{r). A 
shortcoming of this generalization is that it may be not 
unique, but it is the simplest one that provides the points 
mentioned above. 

Equations ( ^6|) and ([l3|) make it possible to express 
e in terms of the pair wave functions and momentum 
distribution. So, a variational procedure can be employed 
to determine these quantities. Perturbing ip(k) and 
and bearing in mind (|4~4|), from Eqs. (|36|) and j43| ) we 
find the following equation 14 1: 



- 2T k i>(k) = U(k) [1 + 2(n k + nrp{k))} , (47) 

where ip(k) and ^ p (k) are the Fourier transforms of the 
scattering waves (|37|). Here U(k) and T k are defined by 
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U(k)= J d 3 r^(r)$(r)exp(-zkr), 
Tk = T k + nt k , 
where t k = U'(k) — U(k) and 



(48) 
(49) 



U'(k)= I d 3 r(^ /2 (r)-^ 2 (r))$(r) 

d 3 g C/(g)(V^ /2 (q)-V- 2 (g)) 



(27T) 



(50) 



Note that t k oc fc 4 for fc -> @, and T fe /T fc -> 1 when 
fc — * 00. So, at low momenta and small densities we 
have ntk <C T k . Owing to this property the difference 
between T k and T k does not play a role when calculat- 
ing the first two terms of the low-density expansions for 
the basic thermodynamic quantities. Therefore, at suffi- 
ciently small values of n Eq. ( |47[ ) is reduced to 

h 2 

— V 2 <p(r) = ip(r)$(r) 
m 

+n f d 3 y<p(yWy)(g tr (\r-y\)-i), (51) 



where gtrif) stands for the truncated pair distribution 
function that is equal to the right-hand side of Ecl(4C]) 
even beyond the weak-coupling regime. Equation (|51) is 
very similar to the Bethe-Goldstone equation pl| . How- 
ever, there is also difference between them. Indeed, g tr 
depends not only on the scattering wave ip{k) but on 
the momentum distribution n k , too. While the Bethe- 
Goldstone equation involves one unknown, the pair wave 
function. 



Equations ( |36[) and ( |44[ ) are accurate to the next-to- 
lcading order in (n—no)/n, then, Eq. ( [47| ) can be accurate 
only to the leading order in (n — no)/n. So, to find ip(k) 
and one should solve Eq. (47) in conjunction with 
Eq. ( [44]) taken in the approximation valid to the leading 
order in the condensate depletion. In other words, one 
has to consider the system of Eqs. ji^ ) and ( f47| ) which 
has the following solution: 




nU(k) 



m = — 2 



2nT k U(k) 
U(k) 



2nT k U(k) 



(52) 
(53) 



Note that as to the scattering states of a pair made of 
one condensed and one uncondensed boson, the goal of 
the present investigation makes it possible not to go into 
details. Below it is only sufficient to limit ourselves to 
the relation ©. 



In the zero-density limit, Eq. (53) is reduced to ^{k) = 
V> (0) (/c) = -t/ (0) (fc)/(27fc), which can be rewritten in the 
form of Eq. (|l5|). So, at sufficiently small densities we 
can express the quantities ijj(k) and n k in terms of the 



vacuum scattering amplitude U^\k) given by Eq. ([l4| ) 
[in the Beliaev's notations U^(k) = /(k,0)]. This is 
totally consistent with the well-known argument of Lan- 
dau Q according to which the thermodynamics of dilute 
quantum gases is determined by the vacuum scattering 
amplitude. Note that the expressions for ?/>(fc) and n k 
derived within the original Bogoliubov model can be ob- 
tained J§|,[l4| from Eqs. fl52j ) and ( |55| ) with replacement 
of T k and U(k) by T k and $(fc), respectively. So, in what 
concerns the expressions for n k and ip(k), the situation 
looks as if we operated with a Bose gas of weakly inter- 
acting quasiparticles with the renormalized kinetic en- 
ergy T k and effective interaction U(r) — cp(r)$>(r). This 
is indeed close to the expectations based on the effective- 
interaction approach of the papers Q . 



VI. SHORT-RANGE BOSON SPATIAL 
CORRELATIONS 

Now, to elaborate on the picture of the short-range 
boson correlations, let us investigate how the correlation 
hole stipulated by the repulsion between bosons at small 
separations changes under the influence of the surround- 
ing bosons. At n — > this hole is completely specified 
by the condensate-condensate pair wave function ip(r), 
which can be found from the definition (E^). Using this 
definition and Eq. (|53|), for the scattering amplitude one 
can find 



[/(fc) = 3(fc) - \ 



d 3 g $(|k-q|)E/(g) 
( 2 ^) 3 . Af 2 + 2nf q U(q) ' 



(54) 



which is the in-medium Lippmann-Schwinger equation. 
Let us rewrite Eq. (|54|) in the form 



U(k) 



m - - 



1 f d 3 q $(\k-q\)U(q) 



(2tt) 3 



T„ 



-I, 



where for / we have 



1=1 1*3- 

2 J (2tt)3 



$(|k-q|)!7(g) $(|k-q|)tT(g) 



2nT q U{q) 

Performing the "scaling" substitution 



qV2mn/fi 



(55) 



in the integral and, then, taking the zero-density limit in 
the integrand, for n — ► we find |M[ 



/ = -a$(fc), 

From Eqs. |34|) and f j5C 

U(k) - C/ (0) (fc) = a$(fc) 
d 3 g $(|k- 



-t/ 3 / 2 (0). 



ir 2 h 3 

it now follows that 



(56) 



(27T) 



2T, 



^[U{q)-U^{q)l (57) 
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where U^Qe) obeys Eq. ( |54| ) with n = 0, i.e. the stan- 
dard Lippmann-Schwinger equation. Introducing the 
quantity = -\U(q) - U^(q)]/(2T q ), for its Fourier 
transform £(r) we find the equation that is nothing else 
but the Schrodinger equation ([l^) with ip^°\r) replaced 
by a + £(?"). As £(r) — > when r — + oo, we can conclude 
that £(r) = aip(°> (r). Hence, for rnOwe get 



V V 71 " 



(58) 



Here y(k,ri) — > 1 when rt — > 0. At A; = the de- 
rived result for, say, the in-medium scattering amplitude 
U(k) coincides with the expansion in na? for the effec- 
tive potential found within the effective-interaction ap- 



proach at the zero temperature (see, e.g., the review 22 



Eq. (4.27)). This shows once more that there are actual 
parallels between our model and approach of Ref. ||. 
However, these parallels are accompanied by significant 
differences. First, the in-medium Lippmann-Schwinger 
equation ( |5^ ) is not a variant of the t matrix equation 
which is frequency dependent, contrary to Eq. ((54). Sec- 
ond, Eq. (|54|) has been found beyond any diagram tech- 
nique by means of the variational procedure. One of its 
important consequences is that the pair wave functions 
that "generate" the in- medium scattering amplitudes [see 
Eq. ( fl8| ) and Eq. ([36]) below] in our approach coincide 
with the pair wave functions that make a contribution 



to g(r) [see Eq. (36)]. By contrast, the modified Bogoli- 
ubov model implies the plane waves for y p (r) (p =/= 0) 
in the pair distribution function (B3) (see Sec. tV), while 



one certainly goes beyond the plane-wave approximation 
when calculating t matrix corresponding to a pair of par- 
ticles with nonzero total momenta. 

With the help of Eq. (jssl); at n — > we obtain the 
following in- medium renormalization for ip(°^ (r) at short 
distances: 



<p(r) ~ <p 



CO) 



(r)(l 



(59) 



Equation ( |59|) is indeed a short-range approximation, 
and this can be understood from the fact that, at suf- 
ficiently large momenta k, the main contribution in the 
integral in the right-hand side of Eq. ( |54|) comes from 
the large momenta q. As ]hn q — f00 U(q)/T q = 0, then 
for k — > oo Eq. ( |54| ) is also (like at n — > 0) reduced to 
the two-body Lippmann-Schwinger equation, and, thus, 
tp(r) obeys the Schrodinger equation ( |15| ) at small boson 
separations. I n ord er to explain the origin of the fac- 
tor C = 1 + 8V no? I ^ in Eq. (|59|), we remind that the 
boundary conditions (]38"|), involved implicitly in the in- 
medium Lippmann-Schwinger equation (Q), are valid for 
r > ro, where rg = n -1 / 3 is the mean distance between 
two particles in the gas. The wave function <p(r), related 
to a couple of bosons in the Bose-Einstein condensate, 
does obey Eq. ( |l5| ) at r — ► 0, but it differs from (p^(r), 
subjected to the boundary condition (|l6|), by the factor 



that is determined by the in-medium effects. Thus, the 
range of validity of Eq. ([59]) is restricted by the region 
r < ro; in the vicinity of ro the behaviour of >f(r) is 
rather complicated; and for r ^> ro the function ip(r) is 
subjected to the following asymptotics: 



ip(r) — ip(r) — 1 ~ — 



1 



2tt 3 / 2 



a 1 



(60) 



The latter evaluation can be obtained from Eq. 
which yields — —\/ira/n/k for k — > 0, and at 

r — > oo we arrive at Eq. (pw, contrary to the two-body 
problem, which implies ip^>(r) ~ —a/r. Thus, the un- 
usual "overscreening" takes place for the wave function 
ip(r) owing to in-medium effects. 

Now we are able to calculate the pair distribution func- 
tion for a dilute Bose gas from Eq. (36). By means of 
the substitution (|55|) in the integral in Eq. ( |36| ) one can 
rewrite the pair distribution function for n — * in the 
form 



9(r) 



1 



n 



(61) 



where the relation ( [iq ) is implied. This is short-range ap- 
proximation for g(r), for the oscillating function cos(kr) 
makes an important contribution to the integral in 
Eq. (^) at large r. To obtain a more concrete infor- 



mation from Eq. (61), one should calculate the conden- 
sate depletion. It can be derived from Eq. ( |52] ) with the 
"scaling" substitution given by Eq. (p5|). This leads to 



n 



^3 

d q n q 
(2tt) 3 n 




(62) 



The result again coincides with that of the effective- 
interaction scheme because the momentum distribution 
( p2| ) is very close to n.& found within the modified Bo- 
goliubov model. Rewriting Eq. (^l]) with the help of 
Eqs. ( ]59"1) and (^), we arrive at the expression ([33]) that 
has earlier been obtained by means of the theorems (^) 
and (p^|). One can see that Eq. ( ^3| ) is also violated at 
the mean distance between particles, ro, since Eq. ( J59| ) 
for ip(r) is broken when r > ro. 

For strongly singular potentials, when ^°)(r = 0) = 0, 
the correct result g(r = 0) = takes place according to 
Eq. d33|). As one can see from (|33|), the correlation hole 
coming from the repulsion of bosons at small particle sep- 
arations gets less marked with an increase of the density 
of the surrounding bosons, which is consistent with usual 
expectations concerning the particle spatial separations. 

Now it is interesting to compare the result given by 
Eq. ( p3|) with the data of the Monte-Carlo calculations for 
the hard-spheres p3j. With the two-body wave function 
( |35l ) for the hard-sphere interaction (|2|), Eq. (|3|) reads 



0, 



r < a 



(1 -a/r) 2 [l + 64\/na 3 /(3 v ^ : )], r > a. 



(63) 
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FIG. 2. The pair distribution function corresponding to 
the hard-sphere potential (|2^) for various values of the gas 
parameter na 3 . Solid lines - the data of Eq. (|63[); □, 
A, v ~ the data of the Monte-Carlo calculations |23| for 
na 3 — 10 -4 , 10~ 2 , 1CP 1 , respectively. The difference becomes 
significant at the mean distance between particles ro, where 
r /a = (na 3 ) -1 ' 3 = 21.5, 4.64, 2.15, respectively. Dotted 
line - Eq. ( fjjf ) for na 3 = 0, one can see that g(r) is essentially 
changed even for small values of the gas parameter. 



The results are displayed in Fig. g for various values of 
the gas parameter na 3 . We remark that, strictly speak- 
ing, only the situation when na 3 = 10~ 4 can be investi- 
gated with Eq. ([33]) because even at na 3 = 10~ 2 we have 
64Vna 3 / w 1.20, so that the next-to-leading order 
in Eq. ( p3| ) makes contribution more essential than that 
of the leading order in na 3 . However, as it is seen, the 
short-range approximation of Eq. (|6^) works well at small 
values of r even at na 3 = 10 — 1 . This agrees with a con- 
clusion of Giorgini et al. ^3] that Eq. ([|), derived to the 
next-to-leading order in na 3 , provides results very close 
to the energy of a dilute Bose gas even at na 3 = 10 _1 . 
Agreement of the data of the Monte-Carlo calculations 
with the results coming from Eq. ( |36| ) can be consid- 
ered as the direct positive test supporting the expression 
$3d(), which is based on the Bogoliubov principle of corre- 
lation weakening in the strongly interacting Bose gas (see 
Ref. [0). _ 

Concluding this section, let us stress one more that 
Eqs. (|5^ ) and (|33|) are correct only for short boson sep- 



arations. To determine the long- and intermediate-range 
behaviour of g[r) with the aim, for example, of finding 
the structural factor S(k) at all values of k, one should 
obtain <p p (r) for all values of p [remember that only the 
limiting relation ( [46"| ) have been utilized]. This inves- 
tigation is directly connected with the problem of the 
relation between the boson momentum distribution and 
pair scattering waves [see Eq. (Q)] and needs additional 
extended considerations being beyond the scope of this 
paper. The same concerns the spectrum of elementary 
excitations. Now, and here the point is, this problem 
is much more complicated as compared to the Bogoli- 
ubov model (original or modified). Within the Bogoli- 
ubov approach the long-range behaviour of g(r) is only 
governed by the two quantities tp(k) — (a^a-kj/no and 
n k = ( a k a k) [see Eq. ( prOf)]. But now we have the rather 
complicated functional (|36|) of nk, ip(k) and ^ p (k) given 
by Eq. m. 



VII. THERMODYNAMICS OF A DILUTE BOSE 
GAS 

After the necessary preparations, we are able to turn to 
the thermodynamics of a strong-coupling Bose gas. The 
most simple way of doing so is to deal with the chemical 
potential u starting from the well-known relation 



u=— L / dV$(|r-r'|)(i/' t (r')V'(r / )V'W} ) 
JnS J 



(64) 



valid in the presence of the Bose condensate pi] . Here 
^T(r) and ip(r) stand for the Bose field operators. This 
relation follows from the expression for the infinites- 
imal change of the grand canonical potential SQ = 
(S(H — uN)) and the necessary condition of the min- 
imum for fl with respect to the order parameter iVo: 
dCl(No, /i, T)/8Nq = 0, the Hamiltonian depending on 
the number of the condensed particles owing to the sub- 



stitution a, 
to 



a 



/iVo. Equations @ and (|§) lead 



n o U{0) + \/2 



where the quantity 



d 3 g 



:%C/ q/2 (q/2), 



E/p(k) = J d 3 r (f p (r)$(r) exp(-ikr) 



(65) 



(66) 



is introduced. The quantities U(k) and t/ p (k), which 
can be called the in-medium scattering amplitudes, differ 
from /(k, k') introduced by Beliaev. The latter is deter- 
mined by the solutions of the ordinary two-body problem 
with the boundary conditions corresponding to the usual 
plane-waves for r — > 00. As to U(k), it is determined 
from the in-medium Lippmann-Schwinger equation (|54|) . 
In order to obtain the next-to-leading order for u, it is 
sufficient to use the relation \/2U(k) — linip^o U p (k), 
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which follows from Eq. (|4(j). In this way, employing the 
substitution (p5|) in the integral and taking into consid- 
eration Eqs. (fil|), ( |58| ) and (p^), we can rewrite Eq. ( |65| ) 
for n — > in the form 



= n£/(Q) ( 1 
47r ft 2 an 




(67) 



Equation ( p7] ) , together with the thermodynamic relation 
jit = d(ne(n))/dn, yields the Yang-Lee's result (Q). 

It is worth noting that we can also use the direct way of 
calculating e based on Eq. ( 143). However, one could come 
to a wrong conclusion if Eq. (|43|) employed in conjunction 
with Eq. (45). To use this way, one should go beyond 
Eq. ( f45| ) and take into account the density correction to 
that relation, 



nk{n k 




n 



n 



V> 2 (fc), 



(68) 



at n — > 0. Thus, the preliminary 
41 should be abandoned in favour 



following from Eq 
result found in Ref. 
ofEq. (|). 

The correct behaviour of the pair distribution function 
found in the previous section allows for deriving the cor- 
rect results for the kinetic and interaction energies by a 
direct calculation, even beyond the Hcllmann-Feynman 
theorem. For example, substituting Eq. ( |33| ) into Eq. (Q) 
and taking account of Eq. (|l^), one can immediately find 
Eq. ©. 



VIII. CONCLUSION 

In conclusion, we remark that the present paper deals 
with the thermodynamics of a Bose gas of strongly inter- 
acting particles in the leading and next-to-leading orders 
of the expansion in na 3 . The strong-coupling general- 
ization of the Bogoliubov model considered here repro- 
duces the well-known formula of Lee and Yang (||) and, 
contrary to the effective-interaction approach of Ref. [Bj , 
yields correct results for the kinetic and interaction en- 
ergies and short-range spatial correlations of bosons. To 
go further, additional investigations should be fulfilled. 
In particular, it is necessary to solve the problem con- 
cerning the relation that connects the boson momentum 
distribution with the scattering waves. The spectrum 
of the elementary excitations should also be considered 
within the approach of Refs. [pf [L3|Jl4| ] in order clarify to 
what extent it differs from the well-known prediction of 
the effective-interaction approach. We, of course, mean 
the region of intermediate momenta rather than the lin- 
ear phonon part, which should be the same according to 
the thermodynamic prescription. This question can not 
be answered without addressing the long-range spatial 
boson correlations. 
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